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We compute the 0(a 3 8 ) total cross section for the process 
pp — > tth in the Standard Model, at ^/~sh = 2 TeV. The next- 
to-leading order corrections drastically reduce the renormal- 
ization and factorization scale dependence of the Born cross 
section and slightly decrease the total cross section for renor- 
malization and factorization scales between m t and 2m t . 



1. Among the most important goals of present and 
future colliders is the study of the electroweak symmetry 
breaking mechanism and the origin of fermion masses. If 
the introduction of one or more Higgs fields is responsi- 
ble for the breaking of the electroweak symmetry, then 
at least one Higgs boson should be relatively light, and 
certainly in the range of energies of present (Tevatron) 
or future (LHC) hadron colliders. The present lower 
bounds on the Higgs mass have been set by LEP to be 
M h > 113.5 GeV fjj] for the Standard Model (SM) Higgs 
boson (h), and M h o A o > 91 GeV for the light scalar 
(h°) and pscudoscalar (A ) Higgs bosons of the minimal 
supersymmctric standard model (MSSM). At the same 
time, precision fits to SM results indirectly point to the 
existence of a light Higgs boson, M h <212- 236 GeV [| , 
while the MSSM requires the existence of a scalar Higgs 
boson lighter than about 130 GeV Q. Therefore, the 
possibility of a Higgs boson discovery in the mass range 
around 115 — 130 GeV seems increasingly likely. 

In this context, the Tevatron will play a crucial role 
and will have the opportunity to discover a Higgs bo- 
son in the mass range between the experimental lower 
bound and about 180 GeV ||. The dominant Higgs pro- 
duction modes at the Tevatron are gluon-gluon fusion, 
gg — > h, and the associated production with a weak bo- 
son, qq — > Wh, Zh . Because of small event rates and 
large backgrounds, the Higgs search in these channels can 
be problematic, requiring the highest possible luminosity. 
It is therefore necessary to investigate all possible produc- 
tion channels, in the effort to fully exploit the range of 
opportunities offered by the available statistics. 

Recently, attention has been drawn to the possibility 
of detecting a Higgs signal in association with a pair of 
top-antitop quarks, i.e. in pp — ► tth o. This produc- 
tion mode can play a role almost over the entire Higgs 
mass range accessible at the Tevatron. Although it has a 
small event rate, ~ 1 — 5 fb for a SM like Higgs, the signa- 
ture (W + W~bbbb) is quite spectacular. Furthermore, at 



the Tevatron (unlike at the LHC), the signal and back- 
ground for this process have quite different shapes. The 
statistics are too low to allow any direct measurement of 
the top Yukawa couplings, but recent studies M indicate 
that this channel can reduce the luminosity required for 
a Higgs discovery at Run II of the Tevatron by as much 
as 15-20%. 

Up to now the cross section for pp — > tth has been 
known only at tree level. As for any other hadronic pro- 
cess, first order QCD corrections are expected to be im- 
portant and are crucial in order to reduce the dependence 
of the cross section on the renormalization and factor- 
ization scales. In this letter we present the results of 
our calculation of the NLO QCD corrections to the total 
cross section for pp — > tth in the Standard Model, at the 
Tevatron. A detailed review of the calculation will be 
presented elsewhere ||. We find good agreement with 
the analogous results presented in Ref. M. 



2. The total cross section for pp 
be written as: 



tth at O{or 8 ) can 



a(pp — > tth) nlo 



^ / dxidx2^Ff(xi,iJ,)^(x2,n) 



■0-NL0(xi,X2,fl) , 



(1) 



where T\ ' p are the NLO parton distribution functions 
for parton i in a proton/antiproton, defined at a generic 
factorization scale /i/ = fj,, and o NLO is the 0(a s ) par- 
ton level total cross section for incoming partons i and 
j, made of the two channels qq, gg — > tth, and renormal- 
ized at an arbitrary scale /x r which we also take to be 
fx r = /i. At the Tevatron, for pp collisions at hadronic 
center of mass energy ^J~sh — 2 TeV, more than 90% of 
the tree level total cross section comes from qq — > tth, 
summed over all light quark flavors. Therefore, we com- 
pute o-(pp — > tth)NLO by including in o NLO only the 
0(a s ) corrections to qq — » tth. The calculation of 
gg — > tth at 0(a^) is, however, crucial to determine 
o~nlo(pP —> tth) for the LHC, since in pp collisions at 
y/s]i = 14 TeV a large fraction of the total cross section 
comes from the gg — > tth channel. The 0(a^) total cross 
section for the LHC has been estimated within the Ef- 
fective Higgs Approximation in Ref. EQ| . Full results are 
presented in Ref. M and will also appear in Ref. Q . 
We write the O(a^) parton level total cross section as: 
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(ari,*2,/x) + <5CT^ LO (xi,x 2 ,/i) , 



where a s (/i) is the strong coupling constant renormalizcd 
at the arbitrary scale // r = /i, <r^ (mi , X% , n) is the O{oq) 
Born cross section, and 6a z ^ LO (xi,X2, fJ.) consists of the 
0(a 8 ) corrections to the Born cross section. The Born 
cross section a l £ (xi,X2, A*) has a strong /^-dependence, 
which is canceled at NLO by Sa l ^ LO (xi,X2, n), up to 
term of 0(a 4 ). The resulting NLO cross section is there- 
fore much more stable under variations of /x, as will be 
discussed in the following (see also Fig. |3|) . 

5ajf LO (xi, X2, fi) contains both virtual and real correc- 
tions to the lowest order cross section and can be written 
as the sum of two terms: 



^tvlo^I'^m) = / d(PS 3 )M(ij -> tth) 

+ j d{PS 4 )M{ij -> tth + g) 

— Cvirt \ &real j 



(3) 



where M(ij — > tth) and M(ij — » ii/l + 5) are respec- 
tively the matrix elements squared for the 0{or s ) 2^3 
and 2^4 processes averaged over the initial degrees of 
freedom and summed over the final ones, while d(PSs) 
and d(PSi) denote the integration over the correspond- 
ing three/four particle phase space. The first term rep- 
resents the contribution of the virtual corrections, while 
the last one is due to the real gluon emission. 

3. The 0(a s ) virtual corrections to the tree level 
qq — > tth process consist of self-energy, vertex, box, and 
pentagon diagrams. The calculation of the virtual dia- 
grams has been performed using dimensional regulariza- 
tion in d = 4 — 2e dimensions. The diagrams have been 
evaluated using FORM |12| and Maple, and all tensor 
integrals have been reduced to linear combinations of a 
fundamental set of scalar integrals. We have computed 
analytically all scalar integrals which give rise to either 
ultraviolet or infrared singularities, while finite scalar in- 
tegrals have been evaluated using standard packages [[13| . 
Among the integrals which have been computed analyti- 
cally, many scalar box integrals and the scalar pentagon 
integrals are extremely laborious, due to the large num- 
ber of massive particles present in the final state. Box 
and pentagon diagrams are ultraviolet finite, but have 
infrared divergences. These integrals are evaluated using 
the method of Rcf. |UJ and analytic results are presented 
in Ref. f§. 

Self-energy and vertex diagrams contain both infrared 
and ultraviolet divergences. The ultraviolet divergences 
are renormalized by introducing a suitable set of coun- 
terterms. Since the cross section is a renormalization 
group invariant, we only need to renormalize the wave 



function of the external fields, the top quark mass, and 
the coupling constants. We use on-shell subtraction for 
the wave-function renormalization of the external fields. 
We define the top mass counterterm in such a way that 
m t is the pole mass. This counterterm must be used 
twice: once to renormalize the top quark mass, and again 
to renormalize the top quark Yukawa coupling. Finally, 
for a s (fi) we use the MS scheme, modified to decou- 
ple the top quark fll5|. T he first nif light fermions are 
subtracted using the MS scheme, while the divergences 
associated with the heavy quark loop are subtracted at 
zero momentum. 

4. The 0(a s ) corrections to the Born cross section 
due to real gluon emission have been computed using 
a two cut-off implementation of the phase space slicing 
algorithm |16|] . The contributions to qq — > tth + g are 
first divided into a soft and a hard contribution, 



2 i(qq — > tth + g) = a soft + a har 



(4) 



where soft and hard refer to the energy of the radiated 
gluon. This division into hard and soft contributions de- 
pends on an arbitrary soft cut-off, S s , such that the en- 
ergy of the radiated gluon is considered soft if E g < <5 S 2 ■ 
The cut-off S s must be very small, such that terms of or- 
der 5 S can be neglected. Therefore, to evaluate the soft 
contribution, the eikonal approximation to the matrix el- 
ements can be taken and the integral over the soft degrees 
of freedom performed analytically. 

The hard contribution to qq — » tth + g is further di- 
vided into a hard/collinear and a hard/not collinear re- 
gion. The hard/collinear region is defined as the region 

where the energy of the gluon is E g > <5 S 2 an( i the 
gluon is radiated from the initial massless quarks at an 
angle 9i g [i — q,q) such that (1 — cos9 g i) < 5 C , for an 
arbitrary small collinear cut-off S c . The matrix element 
squared in the hard/collinear limit is found using the 
leading pole approximation and the integration over the 
angular degrees of freedom is performed analytically. The 
hard gluon emission from the final massive quarks never 
belongs to the hard/collinear region. The contribution 
from the hard/not collinear region is finite and is com- 
puted numerically. 

a so ft and <Jhard/coii contain IR singularities, which 
are calculated using dimensional regularization, and can- 
cel exactly the analogous singularities from the virtual 
contributions, after absorbing mass singularities in the 
renormalized parton distribution functions. 

Both a so ft and (Jhard depend on the two arbitrary cut- 
offs S s and S c , but their sum, i.e. the physical cross sec- 
tion, is cut-off independent. In Figs. [IL we show the de- 
pendence of o~ rea i on the soft cut-off, 5 Sl for a fixed value 
of the hard/collinear cut-off, S c — 10 -4 . In the upper win- 
dow we illustrate the cancellation of the S s dependence 



between o- soft + a hard / coU and a h 



ard/non—coll i 



while in the 



lower window we show a rea i with the statistical errors 
from the Monte Carlo integration. For 8 S in the range 
10~ 4 — 1CP 2 , a clear plateau is reached and the result is 
independent of S s . We point out that Fig. [j] only shows 
distributions of the form da/d8 s , so the corresponding 
cross sections are in this case twice the values that can 
be read from the plot. Analogously, Fig. shows the in- 
dependence of (Treat on the hard/collinear cut-off, S c . All 
the results presented in the following are obtained using 
<L and S r of order 10 -3 . 
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FIG. 1. Dependence of <J rea i(pp — > tth) on the soft cut-off 
<5 s ,at^si"=2TeV,forM h = 120GeV,M = m t ,and5 c = 10~ 4 . 
The lower scale shows the statistical error on a rea i . 
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FIG. 2. Dependence of a rea i(pp — > ith) on the collinear 
cut-off S c , at v /s^" = 2 TeV, for M h = 120 GeV, fi = m t , and 
5 3 =0.005. The lower scale shows the statistical error on a rea i. 



5. Our numerical results are found using CTEQ4M 
parton distribution functions for the calculation of the 
NLO cross section, and CTEQ4L parton distribution 
functions for the calculation of the lowest order cross sec- 
tion Wm. The NLO (LO) cross section is evaluated using 



the 2 (l)-loop evolution of a s (fi). The top quark mass is 
taken to be m t = 174 GeV and af LO (Af z ) = .116. 

First of all, in Fig. | we show, for M h = 120 GeV, how 
at NLO the dependence on the arbitrary renormalization 
scale fj, is significantly reduced. We notice that only for 
scales fj, greater than 2m t is the NLO result greater than 
the lowest order result. 
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FIG. 3. Dependence of <jlo,nlo(pp — * tth) on the renor- 
malization scale fi, at % /sh" = 2 TeV, for M^ = 120 GeV. 



Fig. [| shows both the LO and the NLO total cross 
section for pp — ► tth at ^Js~h — 2 TeV, for two values 
of the renormalization scale fi — m t and /i = 2m t . Over 
the entire range of Mh accessible at the Tevatron, the 
NLO corrections decrease the rate. For example, for 
Mh = 120 GeV and fj, = m t the NLO total cross section is 
reduced to 4.86±0.03 fb from the lowest order prediction 
of 6.868 ± 0.002 fb. The reduction is much less dramatic 
at jjl= 2m t , as can be seen from both Fig. O and Fig. H. 
The complete NLO result includes the lowest order result 
and the contribution of the two terms given in Eq. (H). 
The error we quote on our values is the statistical error 
on the numerical integration involved in evaluating the 
total cross section. 

The corresponding K-factor, i.e. the ratio of the NLO 
cross section to the LO one 



K 



&NLO 
VLO 



(5) 



is shown in Fig. 0. Given the strong scale dependence of 
the LO cross section, the K-factor also shows a significant 
/i-dependence, while it is almost constant with Mh- For 
scales /i between /i = m t and fi = 2m t , the K-factor varies 
roughly between if = 0.70 and K = 0.95. The reduction 
of the NLO cross section with respect to the Born cross 
section is due to fact that at ■JJjj — 2 TeV the tth final 
state is produced in the threshold region. In this region 
the gluon exchange between the final state quarks gives 
origin to Coulomb singularities that contribute to the 



cross section with terms of order a s /fl, where is the 
velocity of the top/antitop quark in the tt CM frame. 
Since the tth final state is in a color octet configuration, 
these corrections are negative and therefore decreases the 
Born cross section, causing the K-factor to be smaller 
than unity. The same effect was observed in the NLO 
cross section for e + e~ — > tth p8j. In that case, however, 
the tth final state is in a color singlet configuration and 
the threshold corrections are positive. 
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FIG. 5. K factor for pp — > tth as a function of Mh, at 
~sh~2 TeV, for /i= mt and /z = 2mt . 



6. The next-to-leading order QCD corrections to the 
Standard Model process pp — > tth, at y/s]j=2 TeV, re- 
duce the cross section by a factor of 0.7 — 0.95 for renor- 
malization and factorization scales m ( </i< 2m t . The 
NLO result shows a drastically reduced scale dependence 
as compared to the Born result and leads to increased 
confidence in predictions based on these results. 
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